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SYMBOLS 


A x amplitude of fundamental component of pressure-coefficient distribution determined by 

Fourier analysis 

A 2 amplitude of second harmonic of the pressure-coefficient distribution determined by 

Fourier analysis 

a wave amplitude, m (in.) 

P " Poo 

Cp pressure coefficient, -g 

H displacement of variable boundary-layer test fixture wall (0 = flush with tunnel wall), m 

(in.) 

L wavelength, m (in.) 

M Mach number 


p static pressure, N/m 2 (lb/ft 2 ) 

p t total pressure, N/m 2 (lb/ft 2 ) 

dynamic pressure, , N/m 2 (lb/ ft 2 ) 

U velocity, m/s (ft/s) 

x distance along wall in direction of air flow, m (in.) 


y distance from wall to position in boundary layer, m (in.) 

P VM 2 - 1 


5 boundary-layer geometric thickness, distance from wall to point where 

U 

fr= 0.98, m (in.) 


6 M=1.0 

5* 

P 


distance from wall to point in boundary layer where M = 1 .0, m (in.) 

Z* 5 / pU \ 

boundary-layer displacement thickness, I (1 — p jf' J dy, m (in.) 
density, kg/m 3 (slugs/ft 3 ) 


0! phase angle of the fundamental component of the surface pressure coefficient relative to 

the wall displacement (positive, leading), deg 


iii 


0 2 phase angle of the second harmonic of the surface static-pressure coefficient relative to 

the wall displacement (positive, leading), deg 

( free-stream condition 
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EXPERIMENTAL INVESTIGATION OF THE INFLUENCE OF THE TURBULENT 


BOUNDARY LAYER ON THE PRESSURE DISTRIBUTION OVER 

A RIGID TWO-DIMENSIONAL WAVY WALL 

Lado Muhlstein, Jr. 

Ames Research Center 

and 

Richard G. Beranek 

George C. Marshall Space Flight Center 


SUMMARY 


The static-pressure distribution was measured on the surface of a series of five rigid sinusoidal 
wavy wall models at Mach numbers from M = 0.80 to 1.35 and ratios of boundary-layer 
displacement thickness to model wavelength of 0.0025 to 0.13. Effects of wavelength and 
amplitude-to-wavelength ratio were examined. For the conditions of the tests, it is shown that the 
turbulent boundary layer causes a large attenuation of the pressure distribution at subsonic speeds 
and both an attenuation and a large phase shift relative to the inviscid case at supersonic speeds. 
These effects are greatest near M = 1 .0 and decrease rapidly with increasing o r decrea sing Mach 
numbers. At supersonic Mach numbers, the pressure coefficient multiplied by >/M 2 — 1 correlates 
closely with the thickness of the subsonic portion of the boundary layer. Results extrapolated to 
zero boundary-layer thickness are presented and compared with three inviscid theories. 


INTRODUCTION 


It has recently been shown that the thickness of the turbulent boundary layer has a large 
influence on the flutter boundary of flat rectangular panels at low supersonic Mach numbers 
(refs. 1-2). The boundary layer has also been suspected of causing large discrepancies between 
theory and experiment in cylinder flutter studies (ref. 3). In addition, the effects of the thickness of 
the boundary layer on the surface pressure distribution are of crucial interest in the study of 
differential ablation or “crosshatching.” As a result, several attempts have been made to study the 
effects of the boundary layer on rigid and moving wavy surfaces (refs. 4-10, for example). In an 
effort to provide basic information that would be useful in developing and evaluating theoretical 
models, a program was undertaken to obtain detailed experimental information on the effect of the 
thickness of the turbulent boundary layer on the pressure distribution on a rigid, two-dimensional, 
sinusoidal wavy wall. Besides being applicable to the fields of panel flutter, cylinder flutter, 
differential ablation, and boundary-layer flow in general, this information, extrapolated to zero 
boundary-layer thickness, should prove useful for evaluating various inviscid flow theories in the 
transonic speed range. 


Prior to this investigation, only one very limited experimental study (ref. 6) was made of flow 
past a wave-shaped wall at low supersonic Mach numbers, and two studies were made of flow past a 
wavy cylinder (refs. 7 and 8). Both wavy cylinder studies showed evidence of flow separation over 
the waves, apparently because the amplitude-to-wavelength ratio was too large. Because the present 
study concentrates on data in the vicinity of M = 1.0 where flow separation is even more likely, it 
was decided to use very small amplitude-to-wavelength ratios. 

Tests were conducted at Mach numbers from 0.80 to 1.35 fora series of five sinusoidal wavy 
wall models and a flat wall model. Wavelengths of 0.05080, 0.15240, and 0.25400 m (2.0, 6.0, and 
10.0 in.), and amplitude-to-wavelength ratios (a/L) of 0.005, 0.010, and 0.015 were studied. 
Results extrapolated to zero boundary-layer thickness are presented and compared with three 
inviscid theories. 


APPARATUS 


Models 

Five sinusoidal wavy wall models were constructed. Three had wavelengths of 0.15240 m 
(6.00 in.) with amplitude-to-wavelength ratios (a/L) of 0.005, 0.010, and 0.015. The other two had 
wavelengths of 0.05080 m (2.00 in.) and 0.25400 m (10.00 in.) with an a/L = 0.005. A typical 
model is shown in figure 1 ; the pertinent dimensions of all the models are given in figure 2 which 
also shows the distribution of the static-pressure orifices. 

For rigidity, the wavy wall models were constructed from 0.076 m (3 in.) thick aluminum 
alloy plate. The sinusoidal wave patterns were formed by a numerically controlled milling machine 
and the resulting surface was hand polished. Measurement showed that all surface points within any 
wave were within 1.27X10' 5 m (±0.0005 in.) of sinusoidal. 


Boundary-Layer Probe 

The probe for measuring the boundary layer on the wavy wall models and the fixture used for 
varying the boundary-layer thickness are shown in figure 3 as installed in the Ames 2- by 2-Foot 
Transonic Wind Tunnel. The probe was designed to traverse the entire length of the model and from 
the surface of the model to 0.076 m (3 in.) into the tunnel stream. Probe position was controlled 
by a programmable servo system so that repeat measurements could be made at the same point. 
Measured reposition accuracy of the system was ±7.62X10' 5 m (±0.003 in.) axially and 
±3.81X10‘ 5 m (±0.0015 in.) normal to the surface of the model. 

The probe was attached to the wind-tunnel model support system with the probe body 
extending into the wind-tunnel diffuser. The probe was fully retracted from the test section while 
surface static pressure was measured and the struts and outboard motor pod were also removed to 
further reduce interference effects. 

The probe was designed for maximum stiffness for positioning accuracy and for minimum 
interference with measured pressures. Maximum tunnel blockage by the probe was limited to 
1.5 percent. 
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The boundary-layer total pressure was measured with a conventional pitot tube with a tip 
opening 1.27X10' 4 m (0.005 in.) high and 1.52X10' 3 m (0.060 in.) wide. The static- pressure 
sensing probe was a 7.62X10" 4 m (0.030 in.) diameter stainless steel tube with a 4° included angle 
conical tip. Two 2.54X10' 4 m (0.010 in.) diameter orifices were located 15 tube diameters 
downstream of the cone-cylinder junction. The static and total pressure probes were insulated from 
the remainder of the probe system and contact with the model surface was indicated electrically. 

Further details of the probe system can be found in reference 1 1 . 


Wind Tunnel 

The experiments were conducted in the Ames 2- by 2-Foot Transonic Wind Tunnel, which is 
the continuous flow type with porous test section walls surrounded by a plenum chamber (see 
ref. 12). Mach number is continuously variable from 0.60 to 1.40 and dynamic pressure is 
continuously variable from 9,576 to 7 1,820 Newtons/M 2 (200 to 1500 psf). Because of the 
installation of the traversing probe in the wind-tunnel diffuser, maximum Mach number was limited 
to 1.35 during these tests. 


Variable Boundary -Layer Test Fixture 

The fixture for varying the boundary-layer thickness consisted of a splitter plate with a sharp 
leading edge installed in one side wall of the wind tunnel. This splitter plate could be positioned by 
remote control to be flush with the wall or to project up to 0.0254 m (1 in.) into the airstream. 
Dimensions of the test fixture are shown in figure 4. When the splitter plate was flush with the 
tunnel wall, the boundary-layer thickness on the model was maximum. When the splitter plate was 
moved into the airstream, a portion of the tunnel wall boundary layer was directed into the plenum 
chamber surrounding the test section and a thinner boundary layer was established on the splitter 
plate. The boundary -layer thickness in the model test region was variable by this means between the 
approximate limits of 0.00254 and 0.03175 m (0.1 and 1 .25 in.) 


INSTRUMENTATION 


Measurement of Pressures 

Wall static pressures and boundary layer rake total pressures were measured by a scanivalve 
system using conventional strain gage differential pressure transducers. Pressure within the 
boundary layer sensed by the traversing probe system was measured by a single pressure transducer. 
The transducer was conncted to the probe tip by a short length of tubing to stabilize the pressure 
rapidly. 
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TEST PROCEDURE 


Flow Conditions in Test Region 

Typical static-pressure distributions on a flat wall in the model test area are shown in 
reference 1. Also shown are typical boundary layer profiles and variations of boundary layer 
thickness with wall position. 


Data Recording Procedures 

The surface static pressures for all models were recorded with the traversing probe in the fully 
retracted position and the struts and outboard pod removed as previously mentioned. 

Detailed surveys of the boundary layers on the wavy wall models were performed at Mach 
numbers of 1.20 and 1.35 by measuring total-pressure profiles and static-pressure profiles during 
separate tunnel runs. The traversing probe could not be used at Mach numbers below 1 .20 because 
of tunnel choking problems and interference with desired results. 

A flat wall model was tested to determine the reference static pressure distribution for use in 
correcting the wavy wall data. These pressures were measured and recorded in the same manner as 
those on the wavy wall. The boundary layer on the flat wall was measured with three total-pressure 
rakes positioned at the four stations shown in figure 4. A different run was required for measuring 
the boundary layer at each station. 


Data Computation 

The surface static pressures measured on the wavy wall models were reduced to coefficient 
form. To improve accuracy, the wavy wall data used in this report were then corrected by 
subtracting the static-pressure coefficient that existed on the flat wall at the same station under 
identical test conditions. The corrections usually varied between Cp = ±0.015 and never exceeded 
±0.030. The corrected surface static- pressure coefficients for each wave with orifices at 1/16 
wavelength intervals were then reduced to their Fourier components. The results are presented as 
the amplitude and phase angle of the fundamental component and the second harmonic. Data from 
the wave farthest downstream usually are not presented because they contained some effects of the 
traversing probe system mounted in the diffuser. 

The boundary-layer velocity and Mach number profiles over the wavy wall models were 
obtained by first reducing the measured static-pressure and total-pressure profiles to coefficient 
form. These data were then used to determine the Mach number and velocity at various heights 
above the model surface. This procedure eliminated the effects of any slight differences in tunnel 
total pressure that may have existed between the static-pressure and total-pressure measurements. 

The boundary-layer velocity and Mach number profiles over the flat wall model were 
computed from the total-pressure profile and a surface static pressure. 
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The height in the boundary layer at which the local Mach number was equal to 1.0 was 
determined from the boundary layer Mach number profile, which was first normalized by the free 
stream Mach number to eliminate the effects of small variations of Mach number during different 
tunnel runs. 


RESULTS AND DISCUSSION 


Influence of Boundary-Layer Displacement Thickness on the Pressure Distribution 

A typical example of the corrected static-pressure distribution over the models is shown in 
figure 5. No separated flow or strong shock waves are apparent. The waveforms are approximately 
sinusoidal. The deviation from a sinusoidal pressure distribution was greatest near M = 1 .0 and for 
models with large a/L because of the increased amount of mixed subsonic and supersonic flow 
present. 

The amplitudes (A } and A 2 ) and phase angles (0! and 0 2 ) of the pressure distributions for all 
five models are plotted in figures 6 to 10 as a function of the dimensionless boundary-layer 
displacement thickness (5*/L). The quantities A t and A 2 are the peak amplitude of the 
fundamental and second harmonic components of the pressure distribution as determined by 
Fourier analysis. The quantities <p l and 0 2 are the corresponding phase angles relative to the wavy 
wall. All data presented in this report were taken on waves with orifices at 1/16 wave intervals. It is 
clear that the thickness of the turbulent boundary layer has a very pronounced effect. At subsonic 
speeds, the effect is primarily an attenuation in amplitude, whereas at supersonic speeds the effect 
is both an attenuation and a phase shift. 

At subsonic speeds 0 X remains near 180° as predicted by inviscid theory and 0 2 remains 
near 90°. Both are essentially unchanged by changes in boundary -layer thickness. At supersonic 
speeds 0! is near 90° for low values of 5*/L and increases with increasing 5*/L. The increase 
in 0i with 5*/L decreases with increasing Mach number; 0 2 behaves similarly except for a much 
greater rate of increase with 5* /L. 

The quantities A l and A 2 decrease with increasing 5*/L at all Mach numbers, with the 
effects greatest at subsonic and low supersonic speeds, and decreasing rapidly above M= 1.20. 
Figure 9(c) at M = 0.95 shows that for 5*/L = 0.1 , A 1 is 66 percent below its value at 5*/L = 0. 

The boundary-layer thickness used in analyzing the wavy wall data is the thickness on the flat 
wall model at the tunnel station corresponding to the center of the wave being analyzed. When a 
wave center did not coincide with a point where the boundary layer was measured, the 
boundary-layer thickness was obtained by interpolation or extrapolation as necessary. 

The boundary-layer data obtained on the wavy wall models were not used in the analysis 
because of the problem of determining the boundary-layer thickness for flow over a wavy wall. 

Although it is expected that the profile shape as well as the thickness of the boundary layer 
will affect the results, no attempt has been made to analyze these data using some measure of 


5 



1 

profile shape as a parameter. Because of the expected effect of the boundary-layer profile shape, an 
integrated boundary-layer thickness (8*) was used rather than the geometric boundary-layer 
thickness (5). The integrated thickness, it is believed, would reduce the sensitivity of the data to 
changes in profiles. In addition, since at many conditions the pressure coefficients measured are 
small, the inevitable resultant scatter in the data would probably mask any systematic variation 
with boundary-layer profile shape. The curves shown in these figures are a cubic polynomial least 
squares fit to the data presented. The extrapolation of these curves to zero boundary-layer 
thickness is shown. Because of the limited amount of data at low values of 8*/L for some models, 
the extrapolation of the least squares fit curve is somewhat uncertain. The steady component (A 0 ) 
is not presented because it was generally small and was not considered significant. Higher harmonics 
were generally very small and erratic except at M = 0.95 where A 3 was approximately 50 percent 
of A 2 • 

In figure 1 1 , the amplitude and phase of the fundamental component and second harmonic 
for models 3—5 are overplotted as a function of 5*/L. The wavelengths of these models are 
0.15240,0.05080, and 0.25400 m (6, 2, and 10 in.), respectively, with a/L = 0.005. This figure 
indicates that, within the accuracy of the data, the quantity 5*/L is a valid dimensionless 
parameter. The curve shown in figure 1 1 is a cubic polynomial least squares fit to the combined 
data. 


The amplitude and phase angles of the fundamental and second harmonic are presented in 
figure 12 as a function of a/L for various constant values of 5*/L including zero. The amplitude 
of the fundamental component A] for 6*/L = 0 is very linear except near M = 1.0. At subsonic 
speeds, this linearity continues with increasing 5*/L whereas at supersonic speeds, Aj as a 
function of a/L becomes less linear with increasing 8*/L. 

The second harmonic A 2 generally increases at an increasing rate with increasing a/L and 
ranges from 10 to 20 percent of A! except at M = 0.95 where it exceeds 30 percent. It should be 
noted that the percentage of A 2 present generally decreases with increasing 5*/L at subsonic Mach 
numbers and increases with increasing 5*/L at supersonic Mach numbers. This occurs because at 
subsonic Mach numbers the boundary layer introduces lower speed flow which is more linear, 
whereas at supersonic Mach numbers a lower speed flow is less linear for a given value of a/L. 

The linearity of A! alone should not be interpreted as proof of aerodynamic linearity since a 
sizable A 2 component is present at some conditions. 

The quantity <t>i is independent of a/L over the entire Mach number range whereas 0 2 is 
independent at subsonic Mach numbers but shows a small dependence on a/L at supersonic Mach 
numbers. 

The amplitude and phase of the fundamental and second harmonic component of the pressure 
distribution as functions of Mach number for various fixed 5*/L are presented in figure 13 for each 
model. Although some of the curves for A], A 2 , and 0 2 for 5*/L = 0 are erratic, the systematic 
behavior of most components at intermediate values of 5*/L indicate that these results are valid. 
The accuracy of the second harmonic component is expected to be much less than the accuracy of 
the fundamental component, particularly at conditions where the second harmonic becomes very 
small. In summary, the extrapolations to 5*/L = 0 should be used with discretion. 
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Figure 13 shows that boundary -layer effects on A 2 are large near M= 1.0 but decrease 
rapidly with increasing or decreasing Mach number. The effect of boundary layer thickness on 
phase angle is small at subsonic speeds and large over the entire supersonic speed range covered by 
these data. An interesting feature shown in figure 13(d) is that, for values of 5*/L — 0.1, the usual 
large increase in A* near M = 1 .0 is almost eliminated. 

The shaded region on the plots of A t in figure 13 indicates the range of 5*/L over which 
experimental data were obtained and is applicable to all four components (A t , A 2 , 0 t , and 0 2 ) for 
each configuration. 

The results of a Reynolds number study performed at M oo = 0.8 are presented in figure 14. 
Amplitude and phase angle of the fundamental and second harmonic components of the static 
pressure coefficients are presented as a function of 5*/L for Reynolds numbers per meter (foot) of 
7.62X10 5 , 1.066X10 6 , and 1.676X10 6 (2.5X10 6 , 3.5X10 6 , and 5.5X10 6 ). All curves are 
essentially identical at all three Reynolds numbers, indicating that S*/L accounts for all Reynolds 
number effects. All other data presented in this report were obtained at a constant unit Reynolds 
number of 1.066X10 6 (3.5X10 6 ) per meter (foot). 

Figure 15 presents the amplitude and phase angle of the combined data from 
the L= 0.05080, 0.15240, and 0.25400 m (2, 6, and 10 in.) models, a/L = 0.005, as a function 
of <5*/L for the four supersonic Mach numbers. These data are plotted in one figure to show the 
large effect of boundary-layer thickness and Mach number. 

This same information is presented in figure 16 in the form of Cp]3 as a function 
of q/L. The quantity 5^ = j q defined as the distance from the model surface to the point 

in the boundary layer over the flat wall where the local Mach number is equal to 1 . The 
determination of this quantity is discussed in the section entitled Data Computation. The pressure 
coefficient amplitudes A t and A 2 were multiplied by j3 to eliminate Mach number effects as 
predicted by linear theory. It is seen that the 5j^_| q boundary-layer thickness very nearly 
collapses all 4> 1 data into a single curve over the full supersonic Mach number range of these data. 
The differences of these curves are within the expected accuracy of these data. The 
quantity 0 2 also collapses, although not as completely as <f> 1 . However, it should be noted that 
the 0 2 data contain a much larger amount of scatter. Multiplying Ai and A 2 by p approximately 
collapses the pressure coefficient amplitude data. These data conclusively demonstrate that at 
supersonic Mach numbers most of the amplitude attenuation and phase shift are produced by the 
subsonic portion of the boundary layer. 

The relationship between 5* and q as derived from the experimental data at the four 

supersonic Mach numbers is shown in figure 1 7. 

Contours of constant Mach number within the boundary layer and into the free stream are 
shown in figure 18 for the models with a/L = 0.005 and L = 0.05080 and 0.15240 m (2 and 6 in.) 
at M = 1 .20 and L = 0.1524 m (6 in.) at M = 1 .35. For reference purposes, 5, 5*, and 5 q for 
the flat wall are shown at each condition. These contours of constant Mach number were derived 
from boundary-layer measurements made with the traversing probe at 17 stations on 
the L = 0.1 5240 m (6 in.) model and at 9 stations on the L = 0.05080 m (2 in.) model. Wall 
static-pressure measurements made with the probe extended and retracted indicated that 
interference at the probe pressure sensing point was small. 


A comparison of the experimental results extrapolated to 5*/L = 0 with theoretical inviscid 
flow results is shown in figure 1 9 as a function of a/L and in figure 20 as a function of Mach 
number. The analytical results used for comparison are linear theory, the nonlinear transonic flow 
theory of Hosokawa (ref. 13) and a nonlinear method of characteristics solution. The values used 
for the linear theory and the nonlinear transonic flow theory of Hosokawa are taken from 
reference 14. The nonlinear method of characteristics results were computed using the program 
discussed in reference 1 5. 

For the conditions where it is valid (fully supersonic flow) the nonlinear method of 
characteristics results agree well with experiment (fig. 19f) for both the amplitude and phase of the 
fundamental component and the amplitude of the second harmonic. The phase angle of the second 
harmonic is in poor agreement; however the experimental values of <t> 2 are known to be of lower 
accuracy. 


CONCLUSIONS 


Based on analysis of the data presented showing the effects of the thickness of the turbulent 
boundary layer on the pressure distribution over a rigid wavy wall, the following can be concluded: 

1 . The thickness of the turbulent boundary layer has a large effect on the amplitude of the 
pressure distribution at all Mach numbers tested. 

2. The thickness of the turbulent boundary layer has a small effect on phase angle of the 
pressure distribution at subsonic Mach numbers, but has a large effect on this quantity at 
supersonic Mach numbers. 

3. The effect of boundary -layer thickness on both amplitude and phase angle of the pressure 
distribution decreases rapidly with increasing supersonic Mach number. 

4. Over the sup ersonic Mach number range of these tests, the pressure coefficient multiplied 
by -y^M 2 — 1 correlates closely with the thickness of the subsonic portion of the boundary 
layer. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., 94035, April 30, 1971 
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Figure 3 -- Front view of traversing probe and variable boundary-layer test fixture installed in 

the Ames 2-by 2-Foot Transonic Wind Tunnel. 





Tunnel station, meters( inches) 



Boundary- layer rake 
Static-pressure orifice 


Note : All dimensions 
are in m (in.) 


Figure 4.- Layout of variable boundary-layer test fixture. 








Model 1 





(a) Moo = 0.80 

Figure 6.- Fourier components of the pressure coefficient on wavy-wall 
model 1 as a function of the dimensionless boundary-layer displace- 
ment thickness; a/L * 0.015, L =* 0.15240 m (6.00 in.). 
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Figure 6.- Continued. 
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Figure 6.- Continued. 
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(a) Moo = 0.80 

Figure 7*“ Fourier components of the pressure coefficient on wavy-wall 
model 2 as a function of the dimensionless boundary -layer displace- 
ment thickness; a/L = 0.010, L = 0.15240 m (6.00 in.). 
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(a) Moo = 0.80 

Figure 8.- Fourier components of the pressure coefficient on wavy-wall 
model 3 as a function of the dimensionless boundary -layer displace- 
ment thickness; a/L = 0.005, L * 0.15240 m (6.00 in.). 
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(a) Moo = 0.80 

Figure 9.- Fourier components of the pressure coefficient on wavy-vall 
model 4 as a function of the dimensionless boundary -layer displace- 
ment thickness; a/L * 0.005, L ■ 0.05080 m (2.00 in.). 
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Figure 11.- Fourier components of the pressure coefficients on wavy- wall 
models 3, 4, and 5 as a function of the dimensionless boundary-layer 
displacement thickness; a/L “ 0.005, L * 0.15240 m (6.00 in.), 

0.05080 m (2.00 in.) and 0.25400 m (10.00 in.). 
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Figure 14.- Fourier components of the pressure coefficients as a function 
of the dimensionless boundary-layer displacement thickness for various 
unit Reynolds numbers at ML * 0.80. 
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Figure 18 .- Constant Mach number contours within the boundary layer and 
free stream above the wavy- wall models. 
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Figure 20.- Comparison of zero boundary-layer thickness experimental 
results as a function of Mach number with predictions of three 
theories . 
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